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Thermal history of the string universe based on the Brandenberger and Vafa's 
scenario is examined. The analysis thereby provides a theoretical foundation of the 
string universe scenario. Especially the picture of the initial oscillating phase is shown 
to be natural from the thermodynamical point of view. A precise description is also 
given of the transition process from the stringy phase to the radiation dominated 
phase. 
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I. INTRODUCTION 



Since the Green and Schwarz's anomaly cancellation [0 proved the importance of the 
super string theory (SST), detailed studies of it have been done. One unfortunate feature of 
SST is the fact that its typical energy scale is 10^^ GeV which is far beyond our experimental 
access. However this does not necessarily imply that SST allows no experimental test. The 
most important feature of SST is that it unifies the theories of matter and the gravity. This 
implies that SST in principle has an ability to determine the structure of space time in which 
the strings themselves live. If SST is the true theory of the whole universe, it is conceivable 
that SST has left some relics in our universe observable even today. In fact the presently 
observed isotropic, uniform and almost flat universe must have been determined by SST. 
From this point of view the string cosmology has been studied by some authors 

Brandenberger and Vafa [§] proposed an interesting scenario of the string cosmology. The 
starting point of their scenario is the Heterotic string theory in the space of nine dimensional 
torus universe of the Plankian size and a time dimension x R. They argued that this 
small universe was oscillating in some period and eventually three dimensions out of nine 
began to expand resulting in the present large universe. 

In order to get deeper theoretical understanding of their scenario, we perform a detailed 
thermodynamical analysis of it in this paper. Our strategy in this paper is to use the 
microcanonical formalizm to follow an entire thermal history. So far several authors have 
employed the microcanonical formalism to examine the thermodynamical functions of the 
string gas. However the relation of their results to the thermal history of the string universe 
seems to be unclear to us. 

In order to clear up these situations in this paper we give a concrete framework by 
which we can follow the thermal history of the string universe using the thermodynamical 
functions of the microcanonical formalism. According to this framework and based on 
some assumptions such as the local thermal equilibrium and others which we will state 
precisely later, we will determine the thermal history of the nine dimensional torus universe 
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of Brandenberger and Vafa as follows. 

In the initial epoch during which the torus universe is oscillating, very high energy strings 
occasionally emit the zero modes (massless point particles) due to the cosmological expansion 
and sometimes absorb the zero modes due to the cosmological contraction. This process is 
shown to be adiabatic. This adiabaticity, or in other words the reversibility, ensures that 
the oscillation is not damping. Thus our result is quite consistent with the picture of the 
initially oscillating universe which is followed by the three dimensionally expanding epoch. 

After some three dimensional directions start to expand with thre remaining six dimen- 
sions being kept Plankian, the energy of zero modes and the strings having no winding along 
the three dimension (which we call non-winding strings) grows roughly in proportion to the 
expanding volume. The temperature is shown to be fixed to be the Hagedorn temperature in 
this period. This inflation like energy growth is possible because the highest energy strings 
(which we call winding strings) continue to supply the energy by their decay. 

This epoch ends when the high energy strings decay away. At this stage there left are the 
dominant zero modes along with a few non-winding strings. From that time the redshift of 
the zero modes get effective resulting in the falling down of the temperature. The remaining 
non-winding string modes are shown to be quickly decay away because of their high specific 
heat. In this way the string universe is shown to transit to the conventional radiation 
dominant universe. The exposition of this thermal history is the main result of this paper. 

Our plan of the discussion is as follows. In the next section we explain on what setting we 
proceed the discussion in this paper and summarize the approximations and the assumptions 



to use in this paper. In Sec. |ITT| we will give a brief review of the Brandenberger and Vafa's 
scenario to fix the notation. In order to follow the thermal history of the string universe, we 
need to calculate the multi string state density. In the ideal string approximation the multi 
state density is evaluated by the single state density. In Sec.^ we give a detailed discussion 
on the single string state density. Especially we explain the change occurring in the single 
state density induced by the cosmological expansion, which is of importance in discussing 
the thermal history. In Sec.0 we give a remark on the value of the total energy of our 



microcanonical ensemble. This value turns out to be the key parameter which determines 
the thermal history. In Sec.[V| we will provide a framework of the microcanonical formalism 
on which we can follow the thermal history. From Sec. |Vll| to Sec.|X| we will evaluate the 



multi state density from the single state densities. In these calculation a novel technique 
is introduced and used extensively. In Sec.|^ the thermal history of our string universe 
is deduced by gathering all the knowledges obtained in the preceding sections. The last 
section is devoted to some discussions. In the Appendendix we will ascertain the validity of 
the Maxwel-Boltzmann approximation which will be use in this paper. 

II. SETTING 

In this section we are going to make sure what kind of tools and approximations we use in 
this paper. As is well known the thermodynamical treatment needs special care for the string 
theory because of the exponentially growing state density One method to treat such 

system is to extend the temperature to complex number in the canonical formalism, and 
the other is to quit to use the canonical formalism and to use the microcanonical formalism 



Both methods are actually connected through the Laplace transformation. We will take 
the latter treatment in this paper. Our interest in this paper is the fundamental string 
theory not a cosmic string. However so as to clarify our setting it is useful to review what 
is known in the studies of the cosmic string. 

The ensemble of the cosmic or fundamental strings is in general subject to both the 
statistical mechanics and the dynamics of the theory. In the case of the cosmic string theory 
the dynamics is shown to prevail over the statistics [jlO|,|lT|. This is a consequence of the 
following settings. First of all for the cosmic string theory one assumes the Einstein gravity 
with Robert son- Walker metric as a background since the relevant energy scale is not so close 
to the Plankian scale. One describes the string as the Nambu-Goto string in the radiation 
or matter dominated back ground. Then the description of the system simplifies well thanks 
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to the one scale principle |]TO|,|TT|. This principle ensures that sooner or later the system will 



be attracted to a scaling solution irrespective of the initial configuration of the strings. It is 
shown that this behavior of the string ensemble is far from the thermal equilibrium. 

However in our case of the fundamental string theory, it is no longer a natural assumption 
that simple Einstein gravity is applicable because the relevant energy scale is as high as the 
Plank mass. The dynamics of the fundamental strings are poorly understood at present, 
thus we cannot proceed further as in the case of cosmic string. That is why we focus on 
the thermodynamical analysis and follow the thermal history using it in this paper. This 
strategy is essentially the one proposed by Brandenberger and Vafa [Q. 

Below we clarify our approximations used in this paper. We are going to investigate 
actually the properties of the string ideal gas in this paper. Thus we need some assumptions 
to identify our system with the real universe. We summarize them here. 

First we have to assume the validity of the local thermal equilibrium and that the ideal 
gas approximation are reasonably good for the string universe in the period of interest. If 
one of these is not a good one, our system of the string ideal gas is not guaranteed to be a 
good approximation of our universe. 

Next we have to assume that some special roles of (quantum) gravity if any are not 
important in considering the thermal history of our universe. Of course string theory by 
birth is the quantum theory unifying gravity. Therefore, some special effects may well exist 
concerning gravity. But our present knowledge is so poor that we have to assume their 
unimportance. 

Even under this difficult circumstance we do think that it is much more meaningful to 
do something rather than doing nothing. Some foothold may well be found by such trial. 
Our main purpose of this paper is to provide the zeroth approximation of the whole story 
of the string universe. 

Within these approximations the thermodynamical functions have been calculated by 
some authors in several models of SST ||^-|TU],|r2|,|I3],0 . In following the thermal history we 
actually need another assumption to determine the history uniquely. As the last assumption 



we require that the usual mechanism of the redshift works for the massless particles even 
in the Plank time. We call this assumption a normal energy loss. In the initial epoch this 
condition is shown to be equivalent to the equi-entropy condition which is also adopted by 



III. COSMOLOGICAL SCENARIO 

In order to fix the notations we present a brief review of the Brandenberger and Vafa's 
I^J^ scenario in this section. 

They started with the Heterotic string theory [|T^ in the nine dimensional torus x R. 
For this model the single string spectrum reads [|^ 
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i=l 



+ {rmaiMl 



(1) 



where 



rrii = rii = 0, ±1, ±2, ■ • • 



nR,nL =0,1,2, 



1/^2^. 



In these expressions \^nai is a linear size of the torus. In the string theory, Ms which is 
of the order of the Plank mass is the only dimensionful constant. We frequently set Mg to 
unity in the sequel. 

The significant feature of this model is a duality a ^ 1/a which is manifest in the 
spectrum. This symmetry connects the large volume world with the small volume world. 
This is called the target space duality The self dual point of this duality is a = 1 (in 



units of 1/Ms). It is known that in the low energy limit of the closed string theory there 
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emerges Einstein gravity [|l|. However Einstein gravity does not respect this duality [Q. 
This means that the use of Einstein gravity is not legitimate in this realm. 

They considered that Einstein gravity is modified in this realm so as to respect this dual- 
ity. The winding mode is thought to play an essential role in this realm. Their approximate 
estimation showed that the winding mode works so as to slow down the expansion of the 
universe while the momentum mode to slow down the contraction. These effects make the 
universe oscillating around the self dual point for a while 

Their description of how three dimensional universe is born is as follows. Suppose that 
a d dimensional space out of nine gets larger than the average by accident. Then the 
winding modes along these d directions get more massive than the rest so that they tend to 
decay more than the average. Consequently, the number of these modes would be reduced. 
Since the winding modes slow down the expansion, these d directions becomes easier to 
expand than the other directions. Thus the perturbation considered above has an unstable 
nature. They suggested that in this way the d dimensional space becomes much larger in 
size than the Plank length while the remaining 9 — c? dimensions are kept Plankian scale. 
Since the Plankian scale is invisible at low energies, this mechanism effectively reduces the 
dimensionality of space time; a, d + 1 < 10 dimensional large universe arises out of small 

X R. They called this mechanism a decompactification. 

IV. SINGLE STATE DENSITY 

In this section we are going to investigate fundamental properties of the single string 
state density f{s). This provides the theoretical foundation of our discussion of the thermal 
history of the string universe. In fact as we see in the later sections, the thermal history 
is deduced from the functional form of the multi state density and the multi state density 
is calculated from the single state density. In the high energy range the functional form 
of f{e) has already been estimated analytically. We first review this result and explain 
how its volume dependence comes out. Later this volume dependence will prove to have 
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key importance in observing the thermal history of the string universe. Next we will give 
our numerical estimation of f{e) by the direct counting of the single string states in the 
low energy range. This clarifies the explicit number distribution of the strings. Lastly we 
remark that an interesting effect in the single state density is induced as the cosmological 
expansion proceeds. 



A. High energy behavior 

For the general closed super string theory in a compact space which is multiply connected, 
the single state density is written as ||7Hl3|,3 

m = l^e^- (2) 

for large enough energy e. In this expression rj = D/2 with D being the number of non- 
compact dimensions, V denotes a D-dimensional volume and 1 / is a constant called the 
Hagedorn temperature. The constants C and (3h depend on the string model. The above 
form of the single state density is uniquely implied by the single string spectrum (0). From 
(0) we learn that the energy of the string consists of a kinetic part r and an oscillation 
part riji + riL. For the kinetic part we have winding modes airrii in addition to the 
usual momentum modes because it is the closed string theory in the non-simply connected 
manifold. 

Before considering the volume dependence of (§) , let us consider what form of the single 
state density is implied in the usual relativistic particle. Such system does not have the 
oscillation mode and the winding mode, the spectrum (|lD reduces to a simpler form = 
J2i=i {^i/aif' on the ci-dimensional torus. The single state density in this case is proportional 
to the surface area of the elliptic sphere in d-dimension having axes eai, ■ ■ ■ ,ead. Namely 
we get 

/(^)«^(n«0^'=^^'"' (3) 

where \^ is a (i-dimensional volume. 



This represents a simple fact that the single state density is an extensive quantity. One 
of the peculiar phenomena in the string thermodynamics is that / is no longer an extensive 
quantity. In fact the number D concerning in (0) is not the total space dimension but the 
noncompact dimension. This means that / is not extensive. In the extreme case of a totally 
compact space, / is volume independent. 

Let us see how this peculiar behavior comes out. Only the kinetic part reflects the 
structure of the space, so that we concentrate ourselves on the degeneracy of the kinetic 
part. Just like the case of the usual point particles, the degeneracy is obtained as the 
surface area of the elliptic sphere having the axes air, ■ ■ ■ , agr, r/ai, ■ ■ ■ , r/ag, see (|l]). The 

state density is therefore proportional to 

d f J, J, J. \ 

— air X x ■ ■ ■ x agr x — x — x ■ ■ ■ x — . (4) 

de V ai a2 ag/ 

This shows that f{e) is certainly a^ independent. This property is essentially a consequence 
of the cancellation between the momentum mode and the corresponding winding mode. As 
the volume expands the phase space of the former increases while that of the the latter 
decreases. Now let us see what happens if D dimensions are open. This time D momentum 
modes miss their partner to cancel, so that a^ dependence remains. Accordingly we get 
f{e) (X ai . . . ao- This surely explains the peculiar volume dependence shown in (|^). 



B. Numerical analysis in low energy range 

We present here the result of our numerical analysis. In Fig.|l| presented is the plot of 
f{e)e~^"^ in the totally compact case D = 0. In the current situation the kinetic energy 
is discrete by the finiteness of the space, which makes the energy spectrum discrete as seen 
in Fig.|l|. We also plot ef{e)e~^'^'^ in Fig.^. The asymptotic behavior f{e)e~^"^ 1/e is 
clearly seen in Fig. 2. This is the first time that this aysmptotic behavior is shown to set in 
already at ~ lOM^. 

These two figures in fact have a clear physical meaning. It is shown |T^ in the micro- 



canonical formalism that f{s)e and ef{e)e represent the number distribution and 



the energy distribution of the strings respectively. 

As we will recognize later, D = 3 case is relevant to our discussion. The value of 
f{e)e~^^'^/V versus e is presented in Fig.^. The spiky behavior therein represents the 
opening of various modes. The analytic estimation indicates that the quantity tends to 
behave C/e^^"^ for large e (see (|^)). The plot e^^'^ f{e)e~^"^ /V is shown in Fig.^ which 
justifies this asymptotic behavior and tells us where this behavior sets in. 

C. Change in the single state density 

In this subsection we examine what occurs to /(e) when accidentally chosen three di- 
rections are expanding while the remaining dimensions are kept Plankian. Because we will 
restrict ourselves to the case in which three directions are expanding at an equal rate, we 
set ai = a2 = = a and 04 = • ■ ■ = ag = 6 ~ 1 (see (|T])) from now on. 

In the preceding section we saw that the high energy behavior of the single state density 
f{e) is independent of a since D = 0. This is the consequence of the cancellation between 
the momentum mode and the winding mode. However this cancellation becomes incomplete 
at low energies for the following reason. 

As a gets larger the winding mode along the a-direction is getting heavier. Eventually the 
winding mode in that direction becomes too heavy to be excited especially in the low energy 
range. This means the winding modes along the expanding three directions are effectively 
frozen. Then the cancellation between the momentum mode and the winding mode breaks 
down in the low energy region. As a result the single state density behaves as if = 3 
instead of = at low energies. Namely f{e)e~^"^ behaves as CV/e^/'^ at low energies and 
as 1/e in the high energy range, respectively. We denote this energy itlq which separates the 
low and high energy ranges. As a becomes large, the strings with higher energies behave as 
if = 3. Namely the effective D = 3 range extends as the universe expands. 

In fact it can be shown by examining the functional form of the state density that as a 
grows mo grows at the rate mo (a) oc 0?. This is justified by the numerical analysis of AUega 
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et al.'s and ours. This phenomena has been discussed also by other authors, in different 
ways by P.Slomonson et al in and by authors in |1^JT^. 



Before closing this section we summarize the behavior of the single state density / in 
an expanding epoch. Below the first excitation energy e < mi, f describes the zero modes 
which is regarded as the usual point particles. In the range between mi and ttiq the state 
density behaves as if the string gas in open three space dimensions. We call the string in this 
range a non-winding string because such string does not have winding along the expanding 
directions. Note that the non-winding string in general has a winding along the remaining 
six directions. 

Lastly for the energy greater than mo(a), / behaves as the string gas in the totally 
compact space. We call such string as a winding string. We can express the single state 
density using the step function 6 as: 

CV 



+ ^e^«"0(£ - mo(a)) 



(5) 



with 1] = 3/2, where fz denotes the state density of the zero modes. Although our interest 
in this paper is in the 1] = 3/2 case only, we keep rj arbitrary in the subsequent expressions 
for better understanding of the structure of our treatment. We again stress here that / is 
proportional to V below mo (a). 



V. HOW BIG SHOULD THE TOTAL ENERGY BE? 

In this section we will give an important remark on the question how big the total energy 
E of our microcanonical ensemble should be. 

Let us consider the string distribution in the initial epoch. The number distribution of 
strings are displayed in Fig.|l|, see the second paragraph of [IV B| . The first excited state 
opens at mi = -\/8(xMs) corresponding to N = nji + = 2 instead of = 1 since the 
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latter is inconsistent with the level matching condition of the Heterotic string theory 

The quantity E is the total energy of the microcanonical ensemble which we have to 
introduce at the beginning of the discussion. What we can find from Fig.|I] is that if we take 
E of the order of Mg as in the usual dimensional analysis, we have no string modes from 
the beginning since the distribution terminates at e = E. In such situation our universe is 
no longer a model of the string cosmology. 

You may say that we only have to take E as large as we like. However in the cosmology 
with causality we can only have finite region in the thermal equilibrium because the speed of 
light is finite. We are dealing with the equilibrium thermodynamics in this paper. Therefore 
it is implicitly assumed that the spatial region having the energy E must be in the thermal 
equlibrium. Consequently we cannot take E as large as we like. 

Now we define E to be the maximally allowed energy in the thermal equilibrium and 
discuss how big E can be. 

Since our present knowledge of the string theory does not allows us to determine the 
value of E, we are left with two possibilities. The first one is that the region having E 
(defined as above) is smaller than the whole universe, the second one is that the whole torus 
universe is in the thermal equilibrium. 

First we consider the former case. Because the value of E is considered to be the maximal 
energy which a single string can occupy, it appears that E must be large enough in order for 
our universe to be regarded as a model of the string cosmology. But it is not necessarily the 
case. In fact there is a loophole in this argument. A string extending over beyond the causal 
region can have the energy much greater than E. Such acausal fundamental string may well 
be produced if the universe itself is born through a quantum tunneling or something like 
that. However such situation is beyond control of our present technology. We do not and 
cannnot go further in such case in this paper. 

Next we consider the second case when the whole universe is in the thermal equilibrium. 
Now we simply conclude that the value of E must be large for the universe to be full of 
strings. There is no loophole this time. 
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Because the loophole in the former possibility is out of our control we decide to assume 
that E is large enough in this paper. Of course the alternative case that our universe has no 
strings even in the initial epoch is another possibility. However we will not treat this case 
since the purpose of this paper is to explore the possibility of the universe full of strings. 

As we stated before our scenario is based on Brandenberger and Vafa's one in which 
the string universe is supposed to oscillate around the self dual point for the initial period. 
This picture nicely fits to the second case mentioned above, since the oscillation over many 
periods tends to thermalize the whole universe. Even if we had started in the loophole case, 
this oscillation makes the acausal strings causal. Moreover, we point out that this picture 
has another advantage from the cosmological view point. This picture is plausible from 
the view point of the horizon problem. If the whole universe would be thermalized during 
the initial period of the oscillation, we would go through the history with the background 
radiation of the same temperature at any part of the universe. 



VI. FRAMEWORK TO FOLLOW THE THERMAL HISTORY 

In this section we intend to give a concrete framework to follow the thermal history of 
the string universe. Many authors have discussed the behavior of the string gas in the micro 
canonical formalism so far [f^-[l3|,^. The str has been examined in various ways and the 



differences lying between the stringy phase (which is frequently referred to as a high density 
phase) and the low temperature phase (which is referred to as a low density phase) have 
been exposed. 

However very little of the attention have been paied to the transition of these two phases. 
The problem how the stringy universe evolves into the radiation dominated universe is still 
an open problem. In order to treat the transient period we present a concrete framework 
based on the microcanonical formalism. 

The multi state density is written by the single state density under the Maxwell- 
Boltzmann (M-B) approximation as |T2[ 
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CO ^ ^ n / " \ 

W = E-t/ I[de,f{e,)s{E-j:eA. (6) 

n=l j=l V i=l / 

We discuss the validity of the M-B approximation in the Appendix. We recall that the 
single state density / is expressed as a sum of the state densities of the zero modes, the 
non-winding strings and the winding strings. We define the multi state densities associated 
with these regions by imitating (H) as: 

oo oo n / " \ 

u!z{ez) = E ZT / n dejfz{ej)5 - E ' 

n=l i=l \ j=l j 

°° I foo n CVde- ( \ 

^N{eN) = E — / n r,+i^ Q^"''0{ej - mi)e (mo(a) - Sj) (5 Uiv - E h 

n=l j=l V j=l J 

°° 1 r°° de ■ ( " \ 

^w{ew) = E — / n {ej - mo(a)) 6 \sw ■ 

n=l ^- -^0 j^i Bj y J 



(7) 

The similarity between (P) and the Taylor expansion of the exponential function enables 
us to anticipate that Q is expressed as a product form of cu's. Actually we can prove by a 
straightforward calculation that 



n{E) = dsz d ENdew uJz{ez)u}N{eN)(^w{£w) 
Jo 

X 6{E -Ez-En - ew)- 

(8) 

The carets on the top of tu's mean the addition of the delta function, uj{e) = uj{e) + 6{e). 
These carets enable us to express the equation in a simple form as above. The necessity 
of the delta functions is readily understood if we recall that right hand side of (P) counts 
the number of all the composite states of three kinds of substances the zero modes, the 
non-winding strings and the winding strings. For example there are also the states having 
no zero modes, which must be counted. The term 6{Ez)ujj^{E]\f)uw{£w) is responsible for 
these states to be taken into account in the integration. 

Based on the equation we argue as follows. Because the delta function ensuring the 
energy conservation is included in the right hand side of (H), we can perform an integral 
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over Ew to obtain a two dimensional integration over {ez, En)- In many cases of interest the 
integrand has a sharp peak at some single point on the {Ez,en) plane, and the contribution 
from this point dominates the integral. The position of the peak is dependent on a since the 
w's have an implicit dependence on a. We denote the position of the peak as (6^(0), eAr(a)). 

If we recall the fundamental principle of the equal a priori probability, we conclude that 
we find the subsystem in the energies (e^(a), cat (a)) when the size of the universe is a. This 
is because this state has an overwhelming probability. This is nothing but the essence of 
the microcanonical formalism. 

Therefore, once we find the functional form of ez{a) and cat (a), we can follow the thermal 
history of the universe. This is our strategy to determine the thermal history in the micro- 
canonical formalism. To carry out this program we need to calculate the multi string state 
densities. In the next three sections we will evaluate the multi state densities associated 
with non-winding strings, winding strings and the zero modes successively. 



VII. MULTI STATE DENSITY OF THE NON- WINDING STRINGS 

The evaluation of the multi state densities ujn{e) and ujiy{e) have been given by several 
authors ||l0 i[T2| , |T3i 3 ■ The method used there is the Laplace transformation and the saddle 



point approximation. We note that the latter is reliable only for large e. However as we 
will see later that the small e behavior of ujn{£) is necessary in the analysis of the thermal 
history. 

In the following we will employ a completely new method to examine the form of ujn 
and uyy. That is the characterization of 's by a differential-difference equation. We will 
show that uj's are solutions to some linear differential-difference equations, and solve them 
to find the functional form of cu's. 
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A. Evaluation 

First we remark that we can factor out the exponential part of U]^ as ujj^^e) = A{e, v)e^"^ 
with 

°° 1 /-OO " . ( " \ 

n=l j=l \ 3=1 / 

where v = CV. This is possible since the integration is constrained by the delta function. 
If we denote 



An{e, v) = Yl -^0{ej - mi)e (mo(a) - 

I " 
x5 e — E 



(10) 



is rewritten as 



oo 



A{e,v) = Y.-An{e,v) (11) 



n=l 



We change the variables and obtain 



1 r°° dx- / " \ 

^"(^' ^) = ^ J n ^^(^.- - "^i/^)*^ irnoia)/6 - x,) 1 - ^ . (12) 



By operating eds and r/f 9^, on it we have 



72 ft 

edeAn{e,v) = -{iqn + l)A„(£,f) H ^A„_i(£: - mi,v) ^A„_i(£: - mo 

and 

?7wa^v4„(e,w) = —^An-i{e-mQ,v), 
mo' 

(13) 

respectively. Here we made use of the fact that mo (a) oc oc v^^'^ (see Sec. [IV C| ) . 



From these we obtain 



16 



'1 + e de + r]vdy)An{e,v) 



n 



-rjnAn {e,v) -\ ^ A„_i {e-mi,v). 

ml 



Summing up this equation over all n we get the equation 



(14) 



'1 + ede + rivdy)A{e, v) = —A{e -Tni,v). 



mi 



(15) 



A{e,v)- 



This is the equation exactly satisfied by A. 

When e is large compared with mi, we can use an approximation A{e — mi^ v) 
midsA^SjV) to rewrite the equation as 

V V 

(1 + ede H 7,mide + rivd^)Aie, v) = —^Aie, v). (16) 

m( m( 

In the case that string has high energy density e/mi » ( which is the case studied 
in refs. [rU,r2|,|T3[, we can neglect the third term of the left hand side of (^) in comparison 
with the second term. The equation reduces to 



;i + ede + rivd.,)A{e,v) 



V 

mi 



A{e,v) 



(17) 



This can be solved with ease to obtain 



A{e,v) 



V 



+1- 



e'l 



X exp 



V ^ 



+ 



V 

rjml 



[e - mi, 



(18) 



with some analytic function g{x). The last step function means that A{e,v) = for e < mi 
by definition. Comparing this with the rj = solution which will be obtained in the next 
section, we can determine g up to x as 



g(x) = + {r]0{r]) + l)x + 



(19) 



where 0{t]) represents the function of t] vanishing at rj = 0. 
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Consequently for the string gas which has the energy En such that v/rje^j^ ~ 0, the multi 
state density is expressed as 



. , cv fcv ^ ' 

u}n[£, v) = const X — ^exp | 7^ + pne 



rjm{ 



(20) 



This reproduces the known result 

We have examined the high density region above. However, we will realize that for our 
string universe the region with low energy density e/mi « v/rjrnl is relevant. Thus we 
have to evaluate ujn{SiV) in this low energy density. Before carrying it out we estimate the 
position of the peak of A and the height of it which are of importance in the thermal history. 

Setting e = ez{a) in (|1^) this equation is reduced to an ordinary differential equation 



1 + T]vd^)A{eN{a),v) 



rjmi 



jfA{eN{a),v) 



(21) 



since d^A vanishes on that point. This ordinary differential equation is readily solved and 
we get the height of the peak: 



A{eNia),v) 



const 



V 



-exp 



rfnii^ 



(22) 



The logarithm of this is nothing but the entropy of the non-winding strings. Therefore, 
the expression simply tells us that the entropy produced when D dimensional volume out 
of nine expands is proportional to the expanding volume (note that v = CV). 

Next we determine the functional form of eN{a). From a numerical determination of 
A{e,v) which we will give later, we see that v/rjmi^ is much greater than 6^(0) /mi. In 



this case we can neglect the second term instead of the third one in ([T6|). This reduces the 
equation to 



-niide + r]vd^ I A{e, v) 



-A{e,v). 



mi'? / mi'' 

By solving it we can determine the functional form of A[e.,v) in this region as 

1 



(23) 



A{e,v) 



exp 



V 



nil rim{ 



V 

rim I 



(24) 



where h is some function to be determined by a boundary condition. But the exphcit form 
of h is unnecessary for our present purpose. The function A{e, v) is maximized when the 
function h{x) is maximal. Denoting the position of the peak of h{x) as x = c, we can express 
the position of the peak of A, eN^a) as 



Numerically it can be expressed as cat (a) = 4.00a'^ + const. Namely cat (a) increases in 
proportion to the volume for large a. 

We solved (^) for small e, at growing values of a. Fig.^ shows the plot of In v) 
versus e. As the energy increases, new modes start to open. This fact is exhibited in the 
low energy behavior of A{e, v) as the emergence of several peaks. We see that the position 
of the peak moves to higher energy as a grows. In high energy region, the spiky behavior 
seen at low energies is smeared, resulting in a smooth curve. Fig.^ and Fig.^ show the plots 
of (eAr(a)/mi)/a'^ and \n{A{eN,v))/{v/rimi) versus a, respectively. The behaviors in both 
figures are consistent with the above analytic estimate. 

To prepare for the later usage we present here the plots of a microcanonical temperature 
of the non-winding strings. This quantity defined as 



measures the rate of the entropy increase due to the energy increase. We show the plots of 
Pn/Ph versus e in FigJ|. 

The global decreasing behavior means that the specific heat is globally positive. The 
point of 6 on which Pm = Ph is where \nA{e,v) peaks. The local spiky behavior around 
e ~ mi means the thermodynamical instability in this region. 




rim{ 



(25) 




(26) 



VIII. MULTI STATE DENSITIES OF THE WINDING STRINGS 



Next we calculate the multi state density of the winding strings. We only have to make 
an analogue of the previous discussion. Since 77 = in the present case we get 
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{l+ede)A{e) = A{e-mo) (27) 

instead of (|T3p. Using a similar approximation we have di;A{e) = 0. Determining the 
normahzation by 

A{mo) = r —e{mo - ei)5{e - = 1/mo, (28) 

J mo Bi 

we finally obtain 

u;w{e) = —ef^^'eie - mo). (29) 
mo 

This reproduces the known result [Q. One particular feature of this functional form is that 
ujw{^)^~^'^'^ has no peak. This is not the case for the non- winding strings and the zero 
modes. 



IX. MULTI STATE DENSITY OF THE ZERO MODES 

In this section we estimate the multi state density of the zero modes. When we enter 
into the thermal history of the string universe, the knowledge on the position of the peak of 
the function uJz{e,v)e~^"'^ will be necessary. We will determine this in the last part of this 
section. 

The multi state density ujz{e,v) in (^ can be rewritten as 

uj^ { e,v) = e^^ 




with an arbitrary positive (3 parameter noting the delta function constraint. 

In order to estimate the power series above, we note two facts here. First the sum of the 
exponential function receives its dominant contribution from k = x term. In fact 

it can be verified that the expression x'^/kl if seen as a function of k has a sharp peak at 
k = X. 
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Next we can ascertain that the power series of uJz{e) can be regarded as ^x^/k\. In 
order for it to be vahd it is sufficient that the important contribution to the integral of e's 
comes from the diagonal region £1 ~ £2 ~ ■ • • ~ £^fc- As we mentioned before fz behaves as a 
positive power of e, see the Sec. |1V A|. This enables us to apply the well-known inequality, 



^ExEi ■ ■ ■ Sk < i^^j to conclude that n/z(^) maximizes on the diagonal region justifying 
our claim. 

00 

Because x of J^x'^/kl corresponds to xq = / fz{£)e~^^de of (|30|), we only have to focus 



on the k = Xo term. Let us examine for which energy e this term does not vanish. For 
the delta function not to vanish, the argument of it must vanish. So let us see what is the 

typical value of J2^=i^j iii the delta function of (|30D. This is estimated as 

00 

k Jdeefzie)e-^' 

<J2^j>=^< > = 

jdefz{e)e-P^ 


00 



deefz{e)e~ 



(31) 



We define here h and n as 



roo 

h{f3,v)= / efz{e)e-^'de, 
Jo 

POO 

n{(3,v)= / fz{e)e-^'de 
Jo 



(32) 



for later convenience. As a result we realize that u!z{s, v)e~^^ has a sharp peak aX e = h{P, v) 
with the height exp(?T,(/5, t>)) with some width q{P,v). 

Making use of a function G{x) having a peak at a; = with a unit width and a unit 
height we can express ujz{e,v)e~^'^ as 

uzie, v)e-^^ = G (^^1^) exp {n{(3, v)) . (33) 

Especially in the case /3 = will be relevant in the later application. We write it here with 
the explicit numerical coefficient 
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h{pH,v) = 65.0a\ 



(34) 



It is worth noting that this argument does not apply to the evaluation of uj]\f{e,v). As 



power of e. The exponential part is irrelevant since it can be factored out as usual. While 
the fact that the remaning part is a negative power of e implies the situation opposite 
to the above case since 1 / {eie2 ■ ■ ■ Sk) maximizes in the boundary region such as ei ~ e, 
~ £2 ~ ■ ■ ■ ~ ^1 ~ ""^i- The is the essence of what called Frautshi Carlitz picture. 

X. THERMAL HISTORY 

In this section we explicitly follow the thermal history of our system using the multi state 
densities calculated in the previous sections. The history which we are going to describe 
below consists of two distinct epochs which we refer to as epoch I and epoch II, respectively. 

We insert cj's obtained above into ^. Among the constituents of fi, uj^ and uy/ have the 
same exponential dependence e'^^^. As for the zero modes, we can formally factor out the 
same exponential using the previous formula ( p3D with [3 = [3h- 

If we insert these ujz, ojn and uw into (||), the exponential factors can be combined to 
form expi^PnE) thanks to the delta function 6{e — sw — — ^n) and then be picked out 
from the integral as 



we saw in (0) the single state density was a product of the exponential part and the negative 




X {e{ew - mo) + S{ew)) 



X 6{E -Sz-sn - ew)- 



(35) 
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A. Oscillating epoch and epoch I 



As we mentioned before we find the system with energies Ez, En, and Ew are determined 
as the position of the peak of the integrand. The functions G and A prefer that Ez and En 
take their most probable values, respectively. The winding string energy Ew is adjusted to 
meet the requirement by the delta function since the integrand of does not have other 
dependence on Ew 

The function for the zero modes strongly favors 

Sz = ez{a) = h{(3H, v) = / d££/,(£)e-^«^ (36) 

•J 

see (1321). While the favorable value for the non- winding strings is similarly determined as 
£n = gn{cl) = miv/{rirnl) + const, which is derived in (pSj). Accordingly the most probable 
value of Ew is determined as Ew = ^w{ci) = E — 6^(0) — cat (a). 

Because both 6^(0) and cat (a) grows as a? for large a, the ratio of cat (a) to e^(a) ap- 
proaches to a constant. Numerically, however, we see that €^{0) « 62(0) from (^) and 
(^). Namely the zero modes is always dominant over the non- winding strings. In view 
of ez{a) in (|36D , we can find that the zero modes are distributed in the canonical distribu- 
tion with the Hagedorn temperature 1/(3h- Namely the temperature in this period is fixed 
to be the Hagedorn temperature. Therefore, the total energy of the zero modes grows in 
proportion to the volume. The reason why it is possible is that the winding strings contin- 
uously supply the energy by decay. The supplied energy is also given to the non-winding 
strings, resulting in the growth of Cat (a) found at ([25|). The conversion of the energy into 
the non-winding strings and the zero modes continues until the winding strings disappear 
(i.e. ew{cL) = 0). We call the period epoch I before their disappearence. 

Here we give an important comment on the change of the total energy E. Generally in 



cosmology the total energy is not a constant quantity [|T8|. For example the energy density 
of the radiation dominated universe scales as 

ez{a)/a^ oc (37) 
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implying ez{a) oc 1/a. This energy loss is attributed to the redshift of the radiation due 
to the expansion of the universe. If the size of the universe is multiplied by a factor a, the 
wavelength is multiplied by the same factor. The radiation loses its energy by this effect. 
The energy lost is given to the gravitational field. We call it a normal energy loss. 

This energy loss is deduced from Einstein gravity. It is true that in the period in question 
Einstein gravity is not a reliable approximation because of the possible string corrections. 
However even in this period, we consider that this energy loss works for the zero mode sector. 
Then we decide to take the normal energy loss as the assumption based on which we follow 
the thermal history of the universe. It is amazing to observe that this normal energy loss 
nicely fits to our scenario. We will show it below. 

In the pure radiation case the normal energy loss is described in the differential equation 
as dez{a)/de = —ez{a)/a. In order to apply this to our case we have to take the existence 
of the other modes into account. The energy exchange between the other modes and the 
zero modes is allowed while the energy loss is only through the redshift of the zero modes 
—ez{a)/a. Therefore the normal energy loss now means 

-^(e^(a) + eAr(a) + ewia)) = ^E{a) = --e^ia). (38) 
da da a 

Adding this equation to the previously given conditions we can uniquely determine the 
functional form of ez{a), e7v(a), ewio-) and E{a) throughout the thermal history. 

In order to reveal what (|3^) means we examine the entropy change of the system under 



this assumption, the entropy of this system is estimated from (|35D as 



S = pHE{a) + n{pH, v) + In A(eiv(a), f ). (39) 

The last term comes from the non- winding strings. Let us first examine the case without it. 
Then the change of 5* due to the growth of a reads 

d^ = _^^eJa^^ SniP,,v) _ (40) 
da a a 

Surprisingly we can show that it vanishes. This is verified by combining the fact that 
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—dj3gn{(3H,v) = h{PH,v) = e^ia) and oc (see (0)). These two imply 

62(0) = j^n{PH,v) which imphes dS/da = 0. 

Consequently, we have proved that as far as we neglect the entropy of the non-winding 
strings, the assumption of the normal energy loss is equivalent to the equi-entropy. Let us 
consider below what this fact means. 

The generation of non-winding strings is rephrased as the unknotting of the winding 
strings along the expanding direction. It was the necessary condition for our universe to exit 
out of the oscillation epoch and entering into the three dimensionally expanding universe 
due to the Brandenberger and Vafa's instability. It is natural to suppose the universe is 
oscillating around the self-dual point, until the condition for the unknotting along some 
three directions is met. 

On the other hand what we proved above is that the process in such epoch is adiabatic. In 
other words the process is reversible. Sometimes winding strings decay by generating the zero 
modes as the universe expands and sometimes the winding strings absorb the energy from 
the zero modes as the contraction proceeds. These processes can repeat themselves because 
they are reversible processes. This situation is naturally identified with the oscillating epoch 
of Brandenberger and Vafa's scenario. 

Once the unknotting in some three directions proceeds enough, the entropy generation 
occurs as we have shown above. This time we cannot go back because the entropy is 
generated. Namely the universe is destined to the three dimensionally expanding universe. 
This is natural as we naively expect that the birth of the three dimensional universe is a 
irreversible process. Consequently, we have recognized that the result of our thermal analysis 
is perfectly consistent with our cosmological scenario thus far. 

The epoch I ends when all the winding strings decay away, in other words when all the 
winding along the three directions unknot. The point a = Oq when the epoch I ends is 
determined by solving the equation E{a) = 62(0) + cat (a). We can easily determine the 
functional form of E{a) from (^) as 
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E{a) = Eo + -{e,{l) - eM) 



(41) 



where Eq = E{1) is the initially given total energy. Using this, (|25| ) and ( p4D we obtain 



ao 



^o + |e,(l) 
|e.(l) + e;v(l). 



1/3 



3^0 



1/3 



(42) 



B. Epoch II 

We call the period a > ao an epoch II. In the epoch II there is no energy supply from 
the winding strings. The non-winding strings and the zero modes compete for the limited 
amount of total energy E{a) in this time. The three functions E{a), ez{a) and ei^{a) in this 
epoch is uniquely determined by the following three conditions: 

e^(a)+eAr(a) (43) 
/3iv(eiv(a),t^) (44) 
= -ie.(a). (45) 

The functions in the second lines are the microcanonical temperature defined as (3z{e, v) = 
dsliaujz{e,v) and (pG]). The second equation is an equi-temperature condition. As we stated 
before the most probable value of the energy is determined as a meeting point of the com- 
petetion between the zero modes and the non-winding strings. The functions (3^ and (3^ 
measure how strongly the respective modes compete for the limited total energy. 

These equations provides the rule for the competition between the winding strings and 
the zero modes. Now the thermal history in the epoch II can be understood qualitatively 
with the knowledge of Pn, which is given Fig.||. 

As a grows, the temperature of the zero modes decreases due to the redshift. For the equi- 
temperature condition to be met, the non- winding strings must cool by the same amount. 
From Fig.^ we can readily find that the non- winding strings have very large specific heat. 
A very little change of the temperature corresponds to a large energy change. Thence the 



E{a) = 

Pziez{a),v) = 

-E{a) 
da 
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decrease of the temperature signifies the violent loss of its energy. For this reason the string 
modes surviving in the epoch II quickly decay away as the expansion proceeds. After the 
strings die away the universe come into the usual radiation dominated universe with usual 
redshift E{a) — ez{a) — 1/a. 

XI. DISCUSSION 

In this paper we have performed the thermodynamical analysis of the string cosmology 
focusing on the Brandenberger and Vafa's scenario. Our analysis was based on the following 
assumptions. 

• The local thermal equilibrium. 

• Ideal gas approximation. 

• Normal energy loss. 

• No unexpected effect due to the non- Einstein correction. 

As a result our analysis has presented the following thermal history of the string universe. 
In the very initial epoch Brandenberger and Vafa's scenario suggested that the universe 
oscillates around the self-dual point of the target space duality. Our analysis has shown that 
the emission of the zero modes from the strings due to the cosmological expansion and the 
absorption of it due to the contraction in the initial universe are adiabatic processes. Then 
these could be repeated, resulting the oscillatory behavior. These observations clarify the 
nature of the oscillating period from the thermodynamical point of view. 

Once the accidental three dimensional expansion is triggered, non-winding strings are 
produced with the entropy production. In this process the winding strings decay producing 
the non-winding strings and the zero modes. During this period the temperature is fixed to 
be the Hagedorn temperature, realizing the inflation-like situation. 
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This period ends when the winding modes are exhausted. At that time the temperature 
begins to fall due to the assumed redshift of the zero modes. The surviving string modes 
quickly die away by this cooling process and there emerged the usual radiation dominant 
universe of the standard big bang cosmology. 

Brandenberger and Vafa also addressed the problem how our space time dimensionality 
is determined to be four within their scenario. In the rest of this paper we consider this 
interesting problem and we propose an alternative idea to determine the dimensionality 
of space time. We first review Brandenberger and Vafa's discussion shortly and give our 
reconsideration next. 

In their scenario, an expanding universe begins as a result of accidental growth of the 
torus universe along some d directions in nine dimensional space. In order for the accidental 
perturbations to grow, the strings should collide with each other frequently so as to diminish 
the number of their winding modes along (i-directions, because the winding modes slow down 
the expansion. If they intersect they would probably unwind. 

On the other hand the string is a two dimensional entity if seen in space time. They 
argued that in order for two strings to collide with finite probability, d+1 must be smaller 
than or equal to 2 + 2 = 4, otherwise space time is too broad for two world sheet to have an 
intersection. This is their derivation of the relation (i + 1 < 4. They said that the universe 
continues to make trials and errors until they learn that less than or equal to three spatial 
directions only become large. When the universe finished the course, d < ?> dimensional 
large universe would have been born. Thus, while they presented the intriguing idea for 
understanding why c/ < 3, they did not succeed in giving compelling reasons why d should 
be 3. 

Now we reconsider their argument. We can not fully agree with their discussion to 
determine the dimensionality because of the following reasons. First, it is true that in the 
point particle case especially in (p^ theory in i?'^, the correlation functions are known |19[ to 
be represented in terms of a random walks in R^. From this, it is rigorously proven that 
the theory is free if + 1 > 5. However this is not the case for the string theory. We have 
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repeated the same analysis as 0^ theory in the hght cone string field theory [2^]. We have 
found that interaction term prevents us from constructing an analogous representation 
as the 0^ theory. 

Second, it is true that the low energy point particles cannot travel to the compact 
direction since the unit momenta are too heavy in such direction. This enforces the point 
particles effectively confined in the d dimensional space. In this case it is only d directions 
that the point particles can utilize to colide. However, we are now concerned with the 
collision of the strings long enough to wind the torus universe many times. One can imagine 
with ease that these strings need not be confined in the d dimensional space and can move 
in any directions. Then what aspects are relevant to the string case? 

Let us recall what occurs to the point particles in the expanding universe. The expansion 
makes the mean separation of particles larger. If the time scale of expansion time get shorter 
than that of interaction time, the interaction is effectively frozen. 

This consideration gives rise to the following interesting possibility. Let us assume that 
the accidental expansion in the d{< 9) directions is always too fast to keep their mutual 
intersections. If it is true, the expansion in the d directions play a negative role with respect 
to the unwinding of the winding strings. 

To put it another way, the winding strings only can use the rest 9 — d dimension to 
unwind. Hence the greater value of 9 — d is favored from the view point of killing the 
winding modes along (i-directions. This implies the inequality opposite in direction from 
that of Brandenberger and Vafa. One may say that this implies that (i = 1 is preferred. 
However the story is not so simple. 

As we have observed in the expression (p^), the entropy produced in the process of the 
unwinding is proportional to the expanding volume. This means that the universe prefers to 
expand as many dimensions as possible. The larger value of d is preferred from the second 
law of thermodynamics, the entropy increase. This effect will compete the preceding one. 
Our idea is that the number of the expanding space dimension is determined to be three as a 
result of this competetion. This idea is not yet formulated on a rigorous ground at present. 
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However we think it is one of the plausible candidates of a mechanism to fix the space time 
dimension. We are planning to make a numerical simulation of the strings to obtain some 
indications to this idea. 
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In this appendix we examine what kind of correction is added when we include the full 
quantum statistics. Rough estimation of it has been done in |jl3|. We will consider this in 
our framework. As a result we can show that the correction in our case is very small. 

It is known [|l^ that the thermodynamical functions in the canonical formalism and the 
micro canonical formalism is connected through the Laplace transformation. If we denote 
the single canonical partition function for the non- winding string as /Ar(/9), this can be 
explicitly written as 

r d(3exp (f^iP)) e"' = A{e, v)e''"'. (46) 

J —too ^ ' 

within the M-B approximation, where dj3 = d(3/2Tii. 

The multi string state density with all the statistical effects ujn{£) is given by the similar 
integral, but the integrand should be changed to 

Z{(3)= n exp -f{r(3) . (47) 



r=l,r:odd 
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Namely the integrand is a product of the functions lf{rP). The equation (^) readily implies 

dpexp(-fN{rp)]ef^' = -A{€/r,v/ry"'^^. (48) 

Here we reminds the readers of a well known formula that the inverse Laplace transform 
of the product of two functions is equal to the convolution of the inverse Laplace transforms 
of the two functions. Namely we have 

deide2(f)i{ei)(l)2{e2)S{e - £i - 62). 

(49) 

Now the repeated use of this formula leads us to the the multi state density with the full 
quantum statistical corrections: 

un{£, v) = lim /q°° deides ■ ■ ■ dsk 

X a{ei,v)a{e?j,v /?>)■■■ Ci{ek,v /k) 



too 

00 



(50) 



where we have set a{e,v) = A{e,v)e^"'^ . All the summation and product indices in this 
appendix mean to run only odd integers unless otherwise stated. 

Now we are going to examine the size of the corrections to the M-B approximation using 
(0). Because a{e, v) does not vanish only for e > mi, the product of (^) is actually a finite 
product. 

For the range knii < e < {k + l)mi with k being an integer we consider the sequence 
of odd numbers such that 1 < Tq < ri < ■ ■ ■ < and Tq + ri + ■ ■ ■ + r/ = fc. Only Tq has 
the possibility to become unity. All the corrections acquired by ujn{£,v) in this range are 
written as the form 

dsro ■ ■ ■ dsri a{ero, V/tq) ■ ■ ■ a{eri,V/ri) 
31 



X S{e- (ro&ro H ^n^rj)- 

(51) 

If we set v/rjml = w and e/rrii = x, a is rewritten as 

a = — exp\h(x — w) + w + (52) 
mi 

Using it and making a variable change = rriiXj ( Here j runs even and odd integer.) 
enable us to rewrite Le~^"^ as 



Le ^-^^ = mi J dxo ■ ■ ■ dxiexp 



J2{h{xj - w/rj) + w/rj - (3Hmi{rj - 1) 

j=0 



5 x-E 



j=0 

(53) 



This is a typical form the corrections added to A{6,v) = a{6,v)e . 

This function is expressed as an integration of a product of functions of the form 

exp [h{xj — w /vj) +w/rj — Pn^iif^j — ■ (54) 

First we consider the case tq = 1. We already know that this peaks at x^"'^ = w + c with 
the height ^^^exp(w), where c is the point such that h'{c) = 0. 

Next we consider the general case rj ^ 1. In this case the function is strongly damped 
by the exponential suppression exp(— /?/^'mi(rj — l)xj). Actually the position of the peak is 
now located around xj^"^ ~ 1. It can be verified by examining where the derivative of the 
exponent of (^4]) changes its sign from positive to negative. The derivative in question is 
written as h'{xj —w/rj) — /5i^mi(rj — 1). We recall here that we had the information on the 
derivative h'{x) because it is essentially the microcanonical temperature examined before, 
see (^). We know that 

h'{x — w) = c?3; In A = miPiy{e, v) 

(55) 
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From our previous numerical analysis (see Fig.^ we know that Pn/Ph is very close to unity 
except x ~ 1, we observe that the above derivative is negative all the range except for the 
very edge Xj ~ 1. This means that the position of the peak is x^""^ ~ 1, and accordingly its 
height is negligibly small. It is no longer the exponentially large like ~ exp{v/rimi). This 
is always true if (ETl) contains the factor with rj ^ 1. Consequently we conclude that any 
corrections to M-B approximation are very small. 
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FIGURES 

FIG. 1. The plot of the single string state density of the totally compact model normalized as 
f{e)e~^"^. This is predicted to tend to 1/e. 

FIG. 2. The plot of ef{e)e~^^^. This shows the predicted asymptotic behavior. 

FIG. 3. The semi-log plot of f{£)e~^"^/V in the case that the space having D = 3 open 
dimensions. 

FIG. 4. The semi-log plot of e'l'^ j(e)&~^"^ jV . This shows the expected asymptotic behavior. 

FIG. 5. The plots of \a.A{£^v) versus e at growing values of a's. The figure (a) is for a = 1, 

(b) is for, a = 1.4, (c) is for a = 1.6, (d) is for a = 2.4 and (e) is the close view around e ~ mi of 
(d). 

FIG. 6. The plot of the favored energy of the non- winding strings cat (a) normalized by 
versus a. This is consistent with the analytic estimate. 

FIG. 7. The plot of In ^(ejv (a), versus a. This is consistent with the analytic estimate. 

FIG. 8. The plots of the micro canonical temperature of the non- winding strings Pn{£,v) 
versus e at growing values of a's. The figure (a) is for a = 1, (b) is for, a = 1.4, (c) is for a = 1.6, 

(d) is for a = 2.4 and (e) is the close view around e ~ mi of (d). 
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